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Abstract 

A new class of nonparametric nonconforming quadrilateral finite elements is in- 
troduced which has the midpoint continuity and the mean value continuity at the 
interfaces of elements simultaneously as the rectangular DSSY element [TJ. The para- 
metric DSSY element for general quadrilaterals requires five degrees of freedom to 
have an optimal order of convergence [3J, while the new nonparametric DSSY elements 
require only four degrees of freedom. The design of new elements is based on the de- 
composition of a bilinear transform into a simple bilinear map followed by a suitable 
affine map. Numerical results are presented to compare the new elements with the 
parametric DSSY element. 

1 Introduction 

There have been many progresses for nonconforming finite element methods for many 
mechanical problems for last decades. Nonconforming elements have been a favorite choice in 
solving the Navier-Stokes equations [El III [131 [H] 111 a stable manner. Also, the nonconforming 
nature facilitates resolving numerical locking [21 [HI [15] in elasticity problems. The lowest 
order nonconforming elements have been especially popular numerical methods because of 
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its simplicity and stability property [6j [13j EH]- The quadrilateral elements can be more 
appropriate than the triangular element for problems with complicated geometry. In this 
paper a class of nonparametric nonconforming quadrilateral elements with four degrees of 
freedom will be introduced which has optimal order of convergence for genuinely quadrilateral 
meshes. 

We briefly review some progresses for nonconforming rectangular or quadrilateral ele- 
ments. Han introduced firstly a rectangular element which assumes five local degrees of 
freedom (DOFs) [8] in 1984. Then in 1992 Rannacher and Turek introduced the rotated Qi 
nonconforming elements with two types of degrees of freedom [13]: the four edge-midpoint 
values DOFs and the four edge integrals DOFs. Chen [5] also used the first type of DOFs for 
the same rotated Qi element. Douglas, Santos, Sheen and Ye introduced a new nonconform- 
ing finite element, which we call the DSSY element in this paper, for which the two types of 
degrees of freedom are coincident on rectangular (or parallelogram) meshes [TJ. One of the 
key features of this DSSY element is that it fulfills the mean value property on each edge, 
which facilitates optimal convergence analyses. Arnold, Boffi, and Falk provided a theory of 
convergence order in quadrilateral meshes [T]. A modified DSSY element was introduced in 
[3J, which requires an additional DOF in order to retain an optimal convergence order for 
genuinely quadrilateral meshes. It seems impossible to reduce the number of DOFs from five 
to four as long as one considers a parametric DSSY-type element on quadrilateral meshes 
and still wants to preserve optimal convergence. 

The aim of this paper is to attempt to extend the spirit of rectangular DSSY element to 
genuinely quadrilateral meshes keeping the mean value property with four DOFs, shifting 
from the parametric realm to the nonparametric one. Our starting point is based on a clever 
decomposition of a bilinear map into a simple bilinear map followed by an affine map [T2~| . 
This approach induces an intermediate reference quadrilateral, where a four DOF DSSY- 
type element can be defined. Then the affine map will preserve P\ and the mean value 
property on each edge. We remark that the quadrilateral element introduced in [12] is of 
only three DOFs, and a similar element was introduced by Hu and Shi [9], but without any 
modification they cannot be used to solve fluid and solid mechanics in stable fashion. 

The paper is organized as follows. In section 2 we review some specific properties of the 
DSSY element. Then using the decomposition of a bilinear map into a simple bilinear map 
followed by an affine map, we introduce a family of quadrilateral elements on an interme- 
diate reference quadrilateral, which is of four DOFs. Based on this, we define a family of 
nonparametric quadrilateral elements. Section 3 is devoted to numerical experiments. The 
performance of the new nonparametric DSSY elements and the parametric DSSY element 
is compared in terms of computation time where the nonparametric DSSY elements show a 
clear advantage over the parametric one. 



A class of nonparametric DSSY elements 



3 



2 Quadrilateral nonconforming elements 

In this section we will introduce a nonparametric DSSY element of four local degrees of 
freedom. First of all let us review the (parametric) DSSY element in brief. 



2.1 The DSSY element 

Let Q be a simply connected polygonal domain in IR 2 and (Th)h>o be a family of shape 
regular quadrilateral triangulations of Q with max^er^ diam(X) = h. Let us denote by Eh 
the set of all edges in Th- For an element K G l~h we denote four vertices of K by Vj for 
j = 1, 2, 3, 4. Also denote the edge passing through Vj-i and v 3 - by ej and the midpoint of e,- 
by nij for j = 1, 2, 3, 4, (assuming vo := V4,) as in Figure [TJ The linear polynomials Z13 and 
/24 are defined in a way that two line equations /13 = 0, £ 24 = pass through mi, m 3 , and m 2 , 
m 4 , respectively. Consider a reference square K — [— 1, l] 2 . We use the similar notations for 
vertices, edges, midpoints of if as those of K such as Vj, Sj, and m, for j = 1, 2, 3, 4. 

Let K G 7^ be any quadrilateral. Then there exists a bilinear map J-^ : K — > X such 
that Fk{K) = K. Notice that J-^ can be written as follows: 

T K {Sc) = Vi + ^-^(v 2 -Vi) + ^-^(v 4 -Vi) + ^ ^ ^( Vl -v 2 + V3-V 4 ). (2.1) 



Set 



where 
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MCff Y = {l,xi,x 2 ,^,(xi) - ^(x 2 )}, Z = 1,2, 



Then the degrees of freedom for the DSSY element can be chosen as either four mean 
values over edges or four edge-midpoint values, which turn out to be identical. In other 
words, the DSSY elements fulfill the mean value property: 



— j vda = v(rh j ), j = 1,2,3,4, Vv e AfC^ 

\ e j\ hi 



DSSY 



(2.2) 



In order to retain an optimal convergence order for any quadrilateral mesh, the parametric 
DSSY element needs an additional element £i£ 2 , and therefore the modified reference element 
reads 

■^ C R S i SY * = {1,^1,^2,^2,^(^1) - £1(^2)}, I = 1,2, 
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with an additional degree of freedom 

v(5i)xiX2 dxidx2 • 
The DSSY element on K is then defined by 



K 



UC DSSY 



■K 



{v | v = v o F K l , v G A/"C ~^ sy *} if K is a true quadrilateral, 
{v\v = vo F K l ,v G A/"C~^ SY } if K is a rectangle, 



where Tk is defined by (12. ip . The global parametric DSSY element is defined by 

MC\ = {vh G L 2 {Vt) | Vh\K G A/"C^ 5sy for K G Th,Vh is continuous at the midpoint of each e G 
A/"C^ g = {v h G A/"C^ | U/j is zero at the midpoint of e G £/, D <9f2}. 



2.2 A Class of Nonparametric DSSY Elements 

We are interested in reducing the five degrees of freedom DSSY element to four, but 
still retaining the mean value property ( 12. 2p . It seems that there does not exist a four- 
DOFs parametric quadrilateral element which has an an optimal order convergence rate and 
mean value property simultaneously. Therefore, we seek a candidate among nonparametric 
elements. 

2.2.1 A closer look at the DSSY element 

For the sake of simplicity and the geometrical property of a basis function, we shall focus 
on (piixt) - <pi(x 2 ). 

Let us denote (pi(xi) — <pi(x2) by for convenience. In the reference domain K, the 
function ^(x) can be factorized as 

^(x) = --{xi - x 2 )(xi + x 2 ){x\ + x% - -), 

6 5 

from which one can realize that ^(x) is the product of three polynomials whose zero-level 
sets consist of the two diagonals of K and one circle x\ + x\ — | = in K. At this point, 
a natural question is whether for any quadrilateral K we may find a function satisfying the 
mean value properties by using the similar geometrical idea as ^(x). 

Among the parametric nonconforming elements in [7], ^(x) = ip o J 7 ^ 1 is not a quartic 
polynomial in general if A is a genuine quadrilateral, that is, if Tk is not an affine map. In 
most cases it is a non-polynomial function. Thus VK X ) "would not be similarly regarded as 
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the product of zero level set functions of three geometrical objects, such as two lines and a 
circle. This seems to be one of the limits of using parametric elements. We will thus divert 
from using parametric elements and investigate in trying to find a suitable four degrees of 
freedom element. 

2.2.2 Intermediate Spaces 

To design such a suitable element, we first decompose the bilinear map Tk into a com- 
position of a simple bilinear map followed by an affine map [12] . For this, let us write 

where Sk '■ K — > K is a simple bilinear map defined by 

and Ar '■ K — > K is an affine map given by 

. V1+V2 + V3+V4 vi-v 2 -v 3 +v 4 _ v 1 +v 2 -v 3 -v 4 _ 
A K {x) = ^ + xi + x 2 . (2.3) 

Here, K = Sk(K) is a quadrilateral with four vertices 

Vi = Vi + s, V 2 = v 2 - S, V 3 = V 3 + S, V 4 = v 4 - s, 



(vi - v 2 - v 3 + v 4 , vi + v 2 - v 3 - v 4 ) 1 (vi - v 2 + v 3 - v 4 ) (2.4) 
1 ( (v 4 - Vi) x (v 3 - v 2 ) 



and 



(V 3 - Vl) X (V 4 - V 2 ) V( V 2 - V l) X ( V 4 - V 3 

In fact, the relations of three mappings and three domains can be interpreted as follows. 
For given domain K £ Th and the reference domain K, Tk is a unique bilinear map such 
that J-^(vj) = Vj for j — 1,2,3,4. If we choose an appropriate vector "s and transform K 
into K by using Is then we obtain a unique affine map Ak from K to K. Actually the vector 
s exists uniquely and is given by (12. 4p . The intermediate reference domain K is very useful 
when we construct certain type of basis functions that have specific characters in K, since 
K is connected to the physical domain K by an affine map not by a bilinear map. Adapted 
to this spirit, we will construct basis functions in K instead of K. 

Remark 2.1. Notice that K is convex if and only if 

\si\ + |s 2 | < 1, (2.5) 
where the equality holds if and only if K degenerates to a triangle. 
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Figure 1: A bilinear map Tk from if to if, a bilinear map Sk from K to if, and an affine 
map Ak from if to if. 

Our strategy is to use the intermediate reference domain if, where the ansatz is to set a 
quartic polynomial similarly to (12. 3p as follows: 

£(x) = -^(x^ 2 (x)Q(x), (2.6) 

where £j(x),j = 1,2, are linear polynomials and Q(x) a quadratic polynomial. We seek a 
quartic polynomial ju(x) fulfilling the mean fa/we property (12.21) in if. Naturally set ^i(x) 
and ^(x) to be linear polynomials such that ^i(x) = and ^(x) = are equations of lines 
passing through v 1; v 3 , and v 2 , v 4 , respectively. Then they are given (up to multiplicative 
constants) by 

^i(x) = xi - x 2 + s 2 - Si, (2.7a) 
F 2 (x) = xi + x 2 + Si + S 2 . (2.7b) 

Recall the Gauss quadrature formula: 
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which is exact for quartic polynomials. An application of this formula simplifies the mean 
value property ( 12.21) into the form 

fe-i) + A»(g2i) - 2/Z(m,-) = 0, j = 1, • • • ,4, (2.8) 

where 

gi = m x - £(u 2 + s), g 2 = mi +£(u 2 + s), 
g 3 = m 2 + £(ui + s), g 4 = m 2 - £(ui + s), 
Is = m 3 + £(u 2 - s), g 6 = m 3 - £(u 2 - s), 
g 7 = m 4 - £(ui - s), g 8 = m 4 + £(u x -s), 

together with m^-, j = 1, • • • , 4, are the twelve Gauss points on the edges. Here, and in what 
follows, we adopt the notations for the standard unit vectors: Ui = ( ) and u 2 



K J V 1 , 

Notice that equations of lines for edges a,-, j = 1, ■ ■ • ,4, are given in vector notation as 
follows: 

ei(t) = mi+t(u 2 +s), e 2 (t) = m 2 +t(ui+s), e 3 (t) = m 3 +t(u 2 -s), e A (t) = m 4 +t(u 1 -s) 

for t G [—1, 1]. Consider the quartic polynomial (I2.6p restricted to an edge ej(t),t G [—1,1]. 
Since £i£ 2 is the product of two linear polynomials which vanishes at the other two end 
points of each edge, one sees that 

4(g 2 ,-i)4(g 2i -i) = IifajMev) = (i - FM&iM&j), (t = Jfj- ( 2 - 9 ) 

A combination of (12.81) and (12.91) yields that (12.21) holds if and only if the quadratic 
polynomial Q satisfies 

Q(g 2j -i) + Q(g2i) - 5Q(m,) = 0, j = l,---,4. (2.10) 

A standard use of symbolic calculation gives the general solution of (12.101) in the following 
form 



2 \ 2 / 2, ' 



Q(x) = ( xx + |s 2 J +(x 2 + |si) -r* + c 



,~ 2_,,_ 2_ 6_ 
(zi + 7S 2 )(z 2 + -si) + — sis 2 
5 5 25 



(2-11) 



with r = "^yf - s i — s 2 f° r arbitrary constant c£l. Here, we assume that the coefficient 
of X\ is normalized. Notice that r is a positive real if K is convex due to Remark 12.11 
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Define, for each eel, 

~ ~ 5~ ~ ~ ~ ~ ~ 

fi(x 1 ,x 2 ;c) = --£i{x 1 ,x 2 )£ 2 {xi, x 2 )Q(xi, x 2 ), 

where l\ and £ 2 are defined by ( 12. 7p and Q by ( 12. lip depending on c as well as s. 

We are now in a position to define a class of nonparametric nonconforming elements on 
the intermediate quadrilaterals K with four degrees of freedom as follows. 

1. K = S K (K); 

2. -Pjj(c) = Span{l,xi,x 2 ,/I(xi,X2;c)}; 

3. = {four edge-midpoint values of K} = {four mean values over edges of K}. 

By the above construction it is apparent that for any element p G fg-(c) the mean ^a/ne 
property holds: 

prr / p da = p{mj), j = 1,2,3,4. 

Moreover, the above class of intermediate nonparametric elements is unisolvent for most of 
c. 

Theorem 2.2. Assume that c is chosen such that + s| + | + c Sisjj 7^ 0. T/ien £/ie 
intermediate nonparametric element (^K, P^(c), zs unisolvent. 

Proof. In order to show unisolvency of the space Span{l, xi, £2, x 2 , c)} with respect to 
the degrees of freedom f(rhj),j = 1, • • ■ ,4, denote the functions 1, xi, £2, and Jl{xi,x 2 ;'c) 
by 0i, 02, 03, and 4 , respectively and also define A = (a^) e M4 x4 (M) by = 0j(m fc ). A 
symbolic calculation shows that det(A) = 16(s^ + s| + | + c S1S2); from which A is nonsingular 
for any s G M 2 if and only if c is chosen such that sf + s| + | + c S1S2 7^ 0. This completes 
the proof. □ 

For c = 0, the quadratic equation Q(x) = denotes the circle with center — | Tf 2 
and radius r. In this case, (12. lip can be easily derived by a geometric argument as follows. 
Indeed, assume that Q(x) = denotes the circle with center c = (^J an d radius r so that 
Q(x) = (x - c) ■ (x - c) - r 2 . Then ( TZTUj) implies that 

(gy_i - c) • (g 2i _i - c) + (g 2j - c) ■ (g 2i - c) - 5(m j - c) ■ (mj - c) = -3r 2 , j = 1, • • • , 4. 
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Arrange these equations as follows: 

(c-^ 2i _i)-(c-^ 2j ) = r 2 , j = !,■■■ ,4, (2.12) 

where the points V2j-i an d V2j are given between g 2 j-i and m,, and g 2j - and m,, respectively, 
explicitly defined as follows: with r\ = y|, 

mi-?7(u2 + s), rj 2 = m 1 + r](u 2 + s), 

m 2 + ?7(ui + s), t) 4 = in 2 - ?7(ui + s), 

m 3 + ?7(u 2 -s), rj 6 = m 3 - r](u 2 - s), 

m 4 - 77(111 -s), t) 8 = m 4 + ?7(ui - s). 

Geometrically, (12.121) is equivalent to saying that the location of c is such that the four 
inner products of the vectors c — Tfoj-i an d c — V2ji f° r 3 = 1, • • • ,4, are equal. It is 
straightforward from the equations (12 . 1 2[) for j = 1 and j — 3 to see that c\ = — ?7 2 s 2 , 
and similarly from those for j = 2 and j — 4 to see that c 2 = —rfsi. Then r = r follows 
immediately. Thus c and r are identical to the center and radius of the circle represented in 
(12. lip in the case of c = 0. 



m = 
m = 

^5 = 
^7 = 



2.2.3 The global nonparametric quadrilateral nonconforming elements 

Turn to the physical domain K. It is straightforward to define the finite elements from 
K to K by using the affine map Ak, which enable to retain the mean value property and 
unisolvency. A class of nonparametric nonconforming elements on quadrilaterals K with 
four degrees of freedom as follows. 

1. K = 7 K {K)- 

2. MC K = Pk(c) = Span{l,xi,x 2 ,/i(2;i,x 2 ;c)}; 

3. Tl K = {four edge-midpoint values of K} = {four mean values over edges of K}, 

where /i(xi,x 2 ;c) is a quartic polynomial defined by /i(xi,x 2 ;c) = Jl o A]^(xi, x 2 ; c) = 
-|£i(xi, x 2 )^ 2 (xi, x 2 )q(x±, x 2 \ c), with 

fx(x) = J x o ^-i(x), £ 2 (x) = 1 2 o ^(x), g(x; c) = Q o ^(x). 

Notice that /i(x; c) can be interpreted as a product of two linear polynomials and one 
quadratic polynomial such that the straight lines ^i(x) = and ^ 2 (x) = are passing through 
Vi, v 3 and v 2 , v 4 , respectively and g(x;c) = is an ellipse which is determined to satisfy 
the mean value properties for /I(x). 
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We now define the global nonparametric DSSY element spaces as follows 

AfC^ p = {vh G L 2 {VL) I Vh\x G NCk for K G 7^,^, is continuous at the midpoint of each e G 
MC^q = {vh G MC 1 ^ | i>h is zero at the midpoint of each eG^fl <9f2}. 

Remark 2.3. S'mce t/iese new finite element spaces have the orthogonal property as in j^, 
clearly the optimal convergence order is guaranteed for solving second-order elliptic problems. 
Also, these will be used as a stable family of mixed finite elements for the velocity fields, com- 
bined with the piecewise constant element for pressure, in solving the Navier-Stokes equations 
^ Q3 [Tffl . The nonconforming nature will solve elasticity problems without numerical lock- 
ing, either FTl\ [73] /. 

Remark 2.4. In practice, the choice c = is recommended since it minimizes the number 
of computations in applying quadrature rules. 

Remark 2.5. We may construct basis functions by the same idea in higher order polynomial 
spaces than the quartic polynomial space. But the latter space would be the best choice in 
using quadrature rules. 

3 Numerical results 

In this section we perform numerical experiments for a simple elliptic problem: 

—Au = f in O, 

u = o on on, 

on the domain Q = (0, l) 2 . The source function / is given so that the exact solution is 

w(x) = sin(7rxi) sin(7ra;2). 

We consider two kinds of elements: the parametric DSSY element MC p h , and the non- 
parametric DSSY elements A/"C^ P with c = and c = 1. Also two types of quadrilateral 
meshes were employed: uniformly ^-dependent quadrilateral meshes as shown in Figure [2] 
are used and the nonuniform quadrilateral meshes depicted in Figure [31 which are generated 
by random perturbations of uniform square meshes. The uniformly ^-dependent quadrilat- 
erals become rectangles if 8 = 0, while they degenerate into triangles if 9 = 1. 

The tables are organized as follows: the parametric nonconforming elements in Tables 1 
and 4, the nonparametric nonconforming elements with c = in Tables 2 and 5, and those 
with c = 1 in Tables 3 and 6; the uniformly ^-dependent quadrilateral meshes in Tables 1-3 
and the nonuniform quadrilateral meshes in Tables 4-6. 
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Figure 2: A quadrilateral triangulation and a quadrilateral with parameter < 9 < 1. 

We test several different ^'s, but the convergence behaviors were quite similar and thus 
we report only the case of 9 = 0.7. 

Numerical experiments were performed with increasing values of c such as 10, 100, 1000, 
and so on. The larger c values are chosen, the slower convergence was observed. We only 
present numerics for the two nonparametric elements with c = and 1. At this point we 
recommend readers to use c = for its simplicity. 

As observed in the uniform mesh the convergence order is optimal for both elements 
and numbers are almost identical. However, there is a big difference in computing time 
between two meshes. The parametric element has the total degrees of freedom, 3N 2 while 
the nonparametric element has 2iV 2 -degrees of freedom. Computing time depends on the 
linear system solvers and the size of the system. By using the conjugate gradient scheme, 
the nonparametric elements are about 1.6 times faster than the parametric element as shown 
in Table [0 



h 


DOF 


\\u - u h \\ 0>Q 


ratio 


\u - u h \ hn 


ratio 


1/4 


40 


0.5284E-01 




0.8532 




1/8 


176 


0.1556E-01 


1.76 


0.4458 


0.94 


1/16 


736 


0.4184E-02 


1.89 


0.2274 


0.97 


1/32 


3008 


0.1096E-02 


1.93 


0.1147 


0.99 


1/64 


12160 


0.2810E-03 


1.96 


0.5756E-01 


0.99 


1/128 


48896 


0.7117E-04 


1.98 


0.2883E-01 


1.00 


1/256 


196096 


0.1791E-04 


1.99 


0.1443E-01 


1.00 



Table 1: Computational Results for MC p h0 with 9 = 0.7. 
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h 


hop 


p — u h\ |o,n 


Tcitio 


\u — Uh\i t u 


±cit io 


1 IA 


94 










1/8 


112 


0.1568E-01 


1.79 


0.4302 


0.93 


1/16 


480 


0.4145E-02 


1.92 


0.2213 


0.96 


1/32 


1984 


0.1084E-02 


1.93 


0.1124 


0.98 


1/64 


8064 


0.2788E-03 


1.96 


0.5659E-01 


0.99 


1/128 


32512 


0.7077E-04 


1.98 


0.2839E-01 


1.00 


1/256 


130560 


0.1783E-04 


1.99 


0.1422E-01 


1.00 



Table 2: Computational Results for A/"C^ P with 9 = 0.7 when c = 0. 



h 


DOF 


\\u - u h \\ ,n 


ratio 


\u - u h \^n 


ratio 


1/4 


24 


0.5840E-01 




0.8486 




1/8 


112 


0.1655E-01 


1.82 


0.4452 


0.93 


1/16 


480 


0.4229E-02 


1.97 


0.2261 


0.98 


1/32 


1984 


0.1102E-02 


1.94 


0.1145 


0.98 


1/64 


8064 


0.2836E-03 


1.96 


0.5760E-01 


0.99 


1/128 


32512 


0.7212E-04 


1.98 


0.2887E-01 


1.00 


1/256 


130560 


0.1819E-04 


1.99 


0.1446E-01 


1.00 



Table 3: Computational Results for AfC^ p with 9 = 0.7 when c = 1. 
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Figure 3: A perturbed nonuniform quadrilateral triangulation 
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^4Q0F,-01 
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0.8663E-02 


2.01 
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1/16 
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1.98 
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48896 
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1.99 
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1.00 


1/256 


196096 


0.9350E-05 


2.00 


0.1183E-01 


1.00 



Table 4: Computational Results for MC p h on the perturbed nonuniform mesh. 
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DOF 
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ratio 
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24 


0.3594E-01 




0.7363 




1/8 


112 
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1.00 


1/16 


480 


0.2290E-02 


1.94 


0.1873 


0.98 


1/32 


1984 
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1.00 


1/64 


8064 
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1.98 
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0.99 


1/128 


32512 
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1.99 
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1.00 
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Table 5: Computational Results for A/"C^ P on the perturbed nonuniform mesh when c = 0. 
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A class of nonparametric DSSY elements 



h 


DOF 


\\u - u h \\ 0t u 


ratio 


\u - Uh\i,n 


ratio 


1/4 


24 


0.3598E-01 




0.7370 




1/8 


112 


0.8752E-02 


2.04 


0.3684 


1.00 


1/16 


480 


0.2290E-02 


1.93 


0.1874 


0.98 


1/32 
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0.5842E-03 


1.97 


0.9398E-01 


1.00 


1/64 
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1.98 


0.4725E-01 


0.99 
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1.99 
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1.00 
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2.00 


0.1184E-01 
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Table 6: Computational Results for on the perturbed nonuniform mesh when 



h 


9 = 0.3 


9 = 0.5 


9 = 0.7 


Perturbed mesh 


1/4 


0.4655 


0.3817 


0.4993 


0.3148 


1/8 


0.6490 


0.6490 


0.6390 


0.6053 


1/16 


0.5778 


0.6000 


0.6000 


0.5971 


1/32 


0.6022 


0.6011 


0.5916 


0.5952 


1/64 


0.6088 


0.6046 


0.5927 


0.5949 


1/128 


0.6008 


0.6030 


0.5944 


0.5933 


1/256 


0.6063 


0.6132 


0.6027 


0.5988 



Table 7: Ratio of computing time t(A/'C^ p )/t(A/'C^ ). 



